LAPTH-1250/08 

Automatised full one-loop renormalisation of the MSSM 
I: The Higgs sector, the issue of tan /3 and gauge invariance 

66" 

O I N. Baro^'', F. Boudjeina^'' and A. Semenov^^ 

^ ; 1) LAPTH, Universite de Savoie, CNRS, 

^ ■ BP 110, F-74941 Annecy-le-Vieux Cedex, France 

^ ^ I 2) Joint Institute of Nuclear Research, JINR, 141980 Duhna, Russia 

^ ■ 

(-H , Abstract 

f^- We give an extensive description of the renormalisation of the Higgs sector of the minimal 

O ■ supersymmetric model in SloopS. SloopS is an automatised code for the computation of one- 

U I loop processes in the MSSM. In this paper, the first in a series, we study in detail the non 

rG ■ gauge invariance of some definitions of tan/3. We rely on a general non-linear gauge fixing 

constraint to make the gauge parameter dependence of different schemes for tan f3 at one- loop 

»v^ , explicit. In so doing, we update, within these general gauges, an important Ward-Slavnov- 

^ ■ Taylor identity on the mixing between the pseudo-scalar Higgs, yl", and the Z". We then 

QQ ' compare the tan/3 scheme dependence of a few observables. We find that the best tan/3 

\^ , scheme is the one based on the decay A" — > t^t~ because of its gauge invariance, being 

^^ ' unambiguously defined from a physical observable, and because it is numerically stable. The 

nJ ' oft used DR scheme performs almost as well on the last count, but is usually defined from 

1^^. . non-gauge invariant quantities in the Higgs sector. The use of the heavier scalar Higgs mass 

f^ ' in Ueu of tan /3 though related to a physical parameter induces too large radiative corrections 

QO ' in many instances and is therefore not recommended. 
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1 Introduction 

Were it not for the radiative corrections to the Ughtest Higgs mass [1] , the minimal supersymmetric 
model or MSSM would have been a forgotten elegant model a long time ago. Indeed, at tree-level 
the mass of the lightest Higgs is predicted to be less than the mass of the Z^ boson, Af^o . That 
would have been a real pity from a model whose most appealing and foremost motivation was to 
solve the hierarchy problem and make the Higgs more natural, beside providing a very good Dark 
Matter candidate. The renormalisation of the Higgs sector of the MSSM is therefore important. 
It is also important because it provides a link to the other parameters of the Standard Model, 
namely all the masses of the particles. It also encodes another parameter that can be seen to 
describe the relative scale of the two vacuum expectation values needed for each Higgs doublet 
of the SM, often referred to as tan/3 and which permeates all the other sectors of the MSSM: 
the gaugino/higgsino sector and the sfermion sector. Many renormalisation schemes or definitions 
of this parameter arc unsatisfactory, as wc will see, mainly because they lack a direct physical 
interpretation or do not correspond to a physical and gauge independent parameter. 

The aim of this paper is to give an extensive description of the renormalisation of the Higgs sector 
in SloopS at one-loop. SloopS is a fully automated code for the one-loop calculation of any 
cross section or decay in the MSSM at one-loop. Although there have been a few studies of the 
renormalisation of the Higgs sector, (see [H [3] for a recent review of the Higgs in supcrsymmetry) 
some performed even beyond the one-loop approximation especially as concerns the mass of the 
lightest CP-evcn Higgs [4j [5], looking at the problem afresh while keeping the issue of gauge 
invariancc in mind, will prove rewarding. Moreover our motivation in developing SloopS was also 
to have a full one-loop renormalisation of all the sectors of the MSSM in a coherent way and 
therefore the study of the Higgs sector is a first step. We will point at the non gauge invariance 
of some definitions of tan/3. Although this has been known, sec for example [6], and pointed 
out at two-loop in the usual linear gauge [7j, most practitioners have kept the usage of some 
non-gauge invariant definitions of tan/3 because of their simplicity at the technical level being 
based on definitions involving two-point function self-energies. With the automatisation of the 
loop calculations, considerations and definitions of tan/3 based on three-point functions (decays) 
are hardly more involved than those based solely on two-point functions describing self-energies, 
including transitions. 

In the approach adopted within SloopS, we strive for an on-shcU, OS, renormalisation scheme 
in particular for tan/3. We rely on a general non- linear gauge fixing constraint to make the 
gauge parameter dependence of different schemes for tan/3 at one- loop explicit. In so doing we 
rederive and update the Ward-Slavnov- Taylor identity on the A^ Z^ / H^W^ mixing in the non- 
linear gauge. We then compare qualitatively and quantitatively the tan/3 scheme dependence of a 
few observables. A^ is the CP-odd Higgs scalar and H^ arc the charged Higgses. Wc find that the 
best tan/3 scheme is the one based on the decay A^ -^ t^t^ because of its gauge invariance, being 
unambiguously defined from a physical observable, and because it is numerically stable. The oft 
used DR scheme performs almost as well on the last count, but is usually defined from non-gauge 
invariant quantities in the Higgs sector. The use of the heavier CP-even scalar Higgs mass in lieu 
of tan /3 though related to a physical parameter induces too large radiative corrections in many 
instances and is therefore not acceptable. It has been argued that the definitions within the Higgs 
sector may be considered universal compared to a definition involving a particular Higgs decay 
for example. However, as stressed in [8], staying within the confines of the Higgs sector and the 
Higgs potential, one faces the issue that many definitions may be basis dependent, as we will see 
this will translate at one-loop into issues about gauge invariance for these definitions. As concerns 
the application to the corrections to the lightest Higgs mass our one-loop treatment is certainly 
not up-to-date, however our motivation is to stress the gauge dependence issues and compare the 
impact of the scheme dependence for tan /3 for many observables starting for those directly related 
to the properties of the Higgses of the MSSM, before reviewing in our forthcoming studies [23] the 



impact of tan/3 on obscrvablcs in the charginos/ncutralinos as well as the sferniion sectors. We 
feel that this issue is of importance as is a consistent one-loop OS implementation. 

The present paper is organised as follows. In Section 2 we review the Higgs sector of the MSSM 
at tree-level. This may, by now, be considered trivial but it is a necessary step before we embark 
on the renormalisation procedure. We also detail this part in order to show what might qualify as 
a physical basis independent observable. Section 3 presents the non-linear gauge fixing condition 
that we use. This includes 8 gauge fixing parameters which are crucial in studying many issues 
related to gauge invariance that are not easily uncovered when one works within the usual linear 
gauge. Section 4 constitutes the theoretical core of our analyses and deals with renormalisation, 
introducing counterterms for the Lagrangian parameters and the field renormalisation constants. 
We expose our renormalisation conditions and update the Slavnov- Taylor identities involving the 
A° — Z^ and H^ — W^ transitions. Section 5 is devoted to defining tan/3. We consider a few 
schemes. Before turning to applications and numerical results we briefly describe how our auto- 
matic code is set-up in Section 6. In Section 7 a numerical investigation of the scheme dependence 
and gauge dependence of these schemes is studied taking as examples loop corrections to Higgs 
masses, decays of the Higgses to fermions and to gauge and Higgs bosons. Section 8 gives our con- 
clusions. The paper contains two appendices. Appendix A details the derivation of Slavnov- Taylor 
identity for the A° — Z'^ transition. Field renormalisation may be introduced at the level of the 
unphysical fields before rotation to the physical fields is performed. Appendix B relates these field 
renormalisation constants on the Higgs fields to the one we introduce directly after the physical 
fields are defined. This may help in comparing different approaches in the literature. 

To avoid clutter we use some abbreviations for the trigonometric functions. For example for an 
angle 9, cos6 will be abbreviated as cg, etc... so that we will from now on use tfj for tan/3. 



2 The Higgs sector of the MSSM at tree-level 

2.1 The Higgs Potential 

As known, see for instance [3], the MSSM requires two Higgs doublets Hi and H2 with opposite 
hypercharge. The Higgs potential in the MSSM is given by: 



V = ml\Hi\^ + rr4\H2\'^ + mj^iHi A H2 + h.c.) (2.1) 

+ lig' + 9'')m\' - \H2\r + 'y\hIh2\' 

with HiAH2=H^H^eab (ei2 - -£21 = 1, £« = 0) . 



The mass terms are all soft masses even if both m^ and m^ contain the SUSY preserving |^p 
term which originates from the F-terms. g,g' are, respectiveley, the SU{2)w and U{1)y gauge 
couplings. Decomposing each Higgs doublet field Hi 2 in terms of its components. 
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the tree-level Higgs potential writes as 
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It is illuminating to express the mass matrices in terms of the tadpoles especially for the pseu- 
doscalar states 
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The requirement that t^i and V2 correspond to the true vacua is a requirement on the vanishing 
of the tadpoles. The tadpoles, by the way, are also a trade-off for m\ and m|. Indeed note that 
expressing everything in terms of T^o , all explicit dependence on m\ and m| has disappeared, 
even in the scalar (CP-even) sector. Note that once the tadpole condition has been imposed 

T^;,, = 0, (2.8) 

we immediately find that in both the charged sector and pseudo-scalar sector, there is a Goldstone 
boson {i.e. a zero mass eigenvalue). This is immediate from the fact that 

dei{NGp)=Q. (2.9) 

The masses of the physical charged Higgs, Mh± and the pscudoscalar Higgs, M^o, arc then just 
set from the invariant obtained from 

Ty{Ngp)^vI + vI=v'', (2.10) 

which is another way of seeing that the combination -y is a proper "observable". Indeed after 
gauging we will find that the masses of the weak gauge bosons are 

Mw± = ^9 V , 

M|c = l{9'+9"y- (2.11) 



Then 



M%, = Tr {Ml«) = -ml^ = m^ + ml , (2.12) 

^ W1W2 

M|± = 7\/2o+M^±. (2.13) 

In Eq. (|2.12p . the first equahty does show an implicit dependence on the ratio of vev (tp), but not 
through mf + ?7i|. The latter must be basis independent, as is the combination m\2/viV2- This is 
to keep in mind. 

It is also interesting to note that for the scalar Higgses, there is a simple sum rule that does not 
involve any ratio of vev's. Indeed, taking the trace of A/?o and call the two physical CP-even 
Higgses /i", with mass M/jO, and i7°, with mass Mho, that would be obtained after rotation, we 
get the sum rule 

M^o + Mfjo = M^o + A/|o . (2.14) 

/i° will denote the lightest CP-even Higgs. Let us as a book-keeping device introduce the angle /3. 
At the moment this is just to help have easy notations: 



Cf3 = — , Sfj ^ — with V ~ Jv\ 4- w| . (2-15) 

The determinant of the scalar Higgses on the other hand gives 

Mlo Mfjo = M'lo A/|o clp . (2.16) 

This shows that if we take M^o, Af^o, M^o as input parameters, we first derive M/jO from Eq. (j2.14p . 
then c^a from Eq. (|2.16p . In general with a set of input parameters M/^o , My^o , M^o , c^a < 1 is 
not guaranteed though. We could of course fix c^a {tf3) and derive Mfjo and My^o which is what 
is usually done. 

The soft SUSY breaking mass parameters ra\ ^ yi can be expressed in terms of the physical quan- 
tities, A/^o, A/^o and cp (as for example derived from Eqs. (|2. 14112. 16)) ): 

m\ = s|A/lo - -C2^A'/|o , (2.17) 
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-2S2/3MI0 , (2.18) 



ml = clMlo + -C2pMlo . (2.19) 



2.2 Basis and rotations 

So far the properties of the physical fields like their masses have been derived without reverting 
to a specific basis. The angle /3 defined in Eq. (|2.15p was just a book-keeping device. Still, to 
go from the fields at the Lagrangian level to the physical fields one needs to perform a rotation. 
This should have no effect on physical observables. This naive observation is important especially 
when we move to one-loop. The rotations we will perform will get rid of field mixing. With the 
tadpole condition set to zero, it is clear that the pseudoscalar and charged scalars eigenstates are 
diagonalised through the same unitary matrix. At tree-level this is defined precisely through the 
same angle /? as in Eq. (|2.15p . 

Ngp = Ui-P) ( ? ) Um , C/(/3) = (^ _'^^ J^ ) , C/t(/3) = Ui-P) . (2.20) 

Call Ty , the tadpole matrix defined as 



T, 



% =1 - ^„ I . (2.21) 



The tadpole is. of course, set to zero. But we will leave this zero there in the notation as we will 
need this when we go to the one- loop counterternis. Then the mass matrices for the CP-even, 
CP-odd and charged scalars write 



M^o = %,+MI„Ngp, (2.22) 


M^i = %, + {Ml„ 


+ M^±)Ngp, (2.23) 


Mlo = %,+ mIoNgp + A/|„[/(/3) (Jo] U{~(i) . (2.24) 


The neutral Higgs is diagonalised through 
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The diagonalisation procedure also produces other, sometimes useful, constraints and relations: 
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= Mlosl_p + Ml,cl^p, (2.26) 
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Note that in the decoupling limit, M^o 3> Mzo, one has in effect decoupled one of the Higgs 
doublet, the other has the properties of the SM Higgs doublet. The decoupling parameter is also 
measured with the parameter cp-a -^ Af|o/Af^o for M^o ^ Mz"- 
Therefore, the mass eigenstates in the Higgs sector are given by 



, - uip)[ ^l) = ( Z !^ V!^ 
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2.3 Counting parameters 

Before we embark on the technicalities of renormalisation and the choice of judicious input pa- 
rameters, it is best to review how we are going to proceed in general and how to make contact 
with the renormalisation of the SM. This will help clarify what are the fundamental parameters 
and which are the physical parameters that can be used for a legitimate renormalisation scheme. 
Moreover since some observables belong to the SM, like the W^, Z'^ masses and the electromag- 
netic coupling constant e which are used as physical input parameters in the OS scheme, isolating 
these three parameters means that their renormalisation will proceed exactly as within the OS 
renormalisation of the SM, see [5] for details. 

In the SM, the fundamental parameters at the Lagrangian level for the gauge sector are g and g' . 
The Higgs potential with the Higgs doublet H. 



v{H) = -fi'^n''n + x{n''nf with 

\{0\nm' = y^O, (2.31) 



furnishes the following: /i^ (the "Higgs mass"), A (the Higgs self coupling) and v (the value of the 
vacuum expectation value). We thus have at Lagrangian level, 5 parameters between the Higgs 
sector and the gauge sector. ^^,X,v are not all independent, v, the vacuum expectation value 
(vev), is defined as the minimum of the potential, this is equivalent to requiring no tadpoles. The 
no tadpole requirement amounts to no terms linear in the scalar Higgs. With the tadpole defined 
as T, we have at tree-level 

T ^ vin"^ - Xv^) ^ . (2.32) 

This requirement is to be carried to any loop level. Out of this constraint, the 5 physical param- 
eters in the OS scheme are e, M^;i/±, M^o, Mj^o, T. At all orders one defines, cw = M^^/i/M^o. 
The latter is not an independent physical parameter. Therefore in the SM a one-to-one map- 
ping between the physical set e, M^r±^ M^o, Mfjo,T and the Lagrangian parameters g, g' ,v, fx, X 
is made. 

In the MSSM, the Higgs sector furnishes m\,m\,m\2 the Higgs doublets soft masses and wi,W2 
the vev of the Higgs doublets . The gauge sector is still governed by the U{1)y and SU{2)w gauge 
couplings .9,.g'. The requirement of no tadpoles from both Higgs doublets, and hence any linear 
combination of them, is also a strong constraint. From these seven parameters in all, the physical 
parameters are usually split between the SM physical On-Shell parameters 

e,Mw±:Mzo, (2.33) 

which are a trade-off for g, g', v'^ ^ vf + wf and the MSSM Higgs parameters 

Mao , T^a , T^o ; "t^" , (2.34) 

which are a trade-off for m^, 7712,771^2,^2 /wi. At tree-level we can set tp = v^jvi but this is, as 
yet, not directly related to an observable. While v can directly be expressed as a physical gauge 
boson mass, the ratio vijvx within the Higgs sector does not have an immediate simple physical 
interpretation. Hence the difficulty with this Lagrangian parameter. One possibility is to trade it 
with the mass of one of the CP-even neutral Higgs through Eq. (|2.16[> . 



3 Non- linear gauge fixing 

In SloopS we have generalised the usual 't Hooft linear gauge condition to a more general non- 
linear gauge that involves, thanks to the extra scalars in the Higgs sector, eight extra parameters 
(a, /?, 5, a), K, p, e, 7). Such gauges within the Standard Model had proved useful and powerful [10|,[9] 
to check the correctness of the calculation. We have also exploited these gauges in the one-loop 
calculation of XiXi ~^ 11 1 Z^l [H] and to corrections to the relic density in [T^]. A 7-parameter 
non-linear gauge- fixing Lagrangian based on the one we introduce here is used in [13j . We can 
extend this non-linear gauge fixing so that the gauge-fixing function involves the sfermions also. 
We refrain, in this paper, from working through this generalisation. 

We will take these gauge fixing terms to be renormalised. In particular the gauge functions involve 
the physical fields. Although this will not make all Green's functions finite, it is enough to make 
all S'-matrix elements finite. The gauge-fixing writes 

£Gf ^ _ 1 f+J7-_ 1 |^2|2_ 1 |^7|2^ (3.1) 
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where 
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F^ = d^Z^^ + ^z (« + e/i°+7^°)G°, 

S2W 

F^ = d^j^. (3.2) 



The ghost Lagrangian CP^ is derived by requiring that the full effective Lagrangian, Cf^ ^ be 
invariant under the BRST transformation. As discussed in [HI, this is a much more appropri- 
ate procedure than the usual Fadeev-Popov approach especially when dealing with the quantum 
symmetries of the generalised non-linear gauges we are using. (^brs'C'^ == therefore implies 



It is very useful to also introduce the auxiliary B-field formulation of the gauge-fixing Lagrangian 
HP^ ^ especially from the perspective of deriving some Ward identities. The gauge fixing can then 
be expressed as 

CpF ^ ^^B+B- + ^|B^|2 + ^IB-'I^ + B-F+ + B+F- + B^ F^ + fi-'i^T . (3.3) 



From the equations of motion for the B-fields we recover the usual C'-^^ together with the condition 



B^ = — -?- (^,; = {£,w ^ ^z T ^.-y}) ■ The anti-ghost, c*, is defined from the gauge fixing functions, we 



write 

feRsc* = B' . (3.4) 

Then the ghost Lagrangian writes as 

^Gh ^ ^(c+s^^^F+ + c~5brsF- + c^^BRsF^ + c^<5brs^^)<5brs>C^'^ . (3.5) 

The Fadeev-Popov prescription is therefore readily recovered, C , but only up to an overall 

function, 6brsC~^^, which is BRST invariant. The latter is set to zero for one-loop calculations. 

Our code SloopS implements this prescription automatically leading to the automatic generation 

of the whole set of Feynman rules for the ghost sector. 

For all applications we set the Feynman parameters ^w,z,'y to one. This allows one to use the 

minimum set of libraries for the tensor reduction. Indeed, Cw,z,7 7^ 1 can generate high rank 

tensor loop functions, that would take much time to reduce to the set of scalar functions. 

It is important to stress, once more, that since we do not seek to have all Green's functions finite 

but only the S-matrix elements, we take the gauge fixing Lagrangian as being renormalised. 

Judicious choices of the the non-linear gauge parameters can lead to simplifications like the van- 
ishing of certain vertices. For example, with (5 = 1. the VF"*"^G~7^ vertex cancels. More examples 
can be found in Appendix A for the vanishing of some ghost couplings to Higgses. 



4 Renormalisation 

Our renormalisation procedure is within the spirit of the on-shell scheme borrowing as much as 
possible from the programme carried strictly within the Standard Model in [S]. For the gauge 
sector and the fermion sector, beside the electromagnetic coupling which we fix from the Thomson 
limit, we take therefore the same set of physical input parameters, namely the masses of the 
W^ and Z^ together with the masses of all the standard model fermions. To define the Higgs 
sector parameters, the set of Eq. p.34p looks most appropriate were it not for the ill defined tp. 
Indeed, the mass of the pseudoscalar Mj^o within the MSSM with CP conservation is a physical 
parameter. As within the Standard Model, we also take the tadpole. For t^ the aim of this paper 
is to review, propose and compare different schemes. Renormalisation of these parameters would 
then lead to finite S-matrix elements. For the mass eigenstates and thus a proper identification of 
the physical particles that appear as external legs in our processes, field renormalisation is needed. 
S-matrix elements obtained from these rescaled Green's functions will lead to external legs with 
unit residue and will avoid mixing. Therefore one also needs wave function renormalisation of the 
fields. Especially for the unphysical sector of the theory, the precise choice of the fields redefinition 
is not essential if one is only interested in S-matrix elements of physical processes. It has to be 
stressed that one can do without this if one is willing to include loop corrections on the external 



legs. In the MSSM and in the Higgs sector in particular mixing effects, especially at one-loop, are a 
nuisance that has introduced some confusion especially in defining tp with the help of wave- function 
renormalisation constants or cquivalently from two-point function describing particle mixing. For 
the Higgs sector one needs to be wary about mixing of the Goldstones with CP-odd Higgs or 
almost cquivalently between the Z^ and the CP-odd Higgs or the W^ and the charged Higgs. 
These two-point functions are related through gauge invariance and impose strong constraints 
on the wave function renormalisation constants. We will derive Ward-Slavnov- Taylor identities 
relating these two-point functions, and hence their associated countcrterms, before imposing any 
ad-hoc condition. 

4.1 Shifts in mass parameters and gauge couplings 

All fields and parameters introduced so far are considered as bare parameters with the exception 
of the gauge fixing Lagrangian which we choose to write in terms of renormalised fields. Care 
should then be exercised when we split the tree- level contributions and the countcrterms. Shifts 
are then introduced for the Lagrangian parameters and the fields with the notation that a bare 
quantity is labeled as Xq. It will split in terms of renormalised quantities X and countcrterms 6X 

90 = 9 + 69, 90 ^9' + 69', (4.1) 

™?o = mf + Smf fori = 1,2, mj^ = mj^ + 5ml2 , (4.2) 
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Vio = Vi — Svi for i = 1,2 hence — - = . (4.3) 

tp Vl V2 

In our approach the angles defining the rotation matrices, P and a in Eq. (j2.30p are defined 
as renormalised quantities. For example the relation between the Goldstone boson/pseudoscalar 
Higgs boson and the fields tpi 2 i^ maintained at all orders. Indeed, 

^M ---^^n implies also (%)-UiP)(^l 



^0 I = Ui(3) "2 impHes also (40 ) - U{P) ( ^J ) ■ (4-4) 



Since in our approach we will always perform a field renormalisation there is no need in inducing 
more shifts from U{a,/3). Therefore U{a,f3) is taken as renormalised. For example, if we perform 
a field renormalisation in the ipi 2 basis 
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(4.6) 

For the field renormalisation we can perform this either at the level of the 1^^, i.e. before any 
rotation on the field in the Lagrangian is made, through Z^o as is done in [1411151 [T5] or in a much 
efficient way directly in the basis G^ J^ since the latter are directly related to our renormalisation 
conditions on the physical fields as we will see later. For instance, there is no need for Zqoqo in 
our approach since we will not be dealing with Goldstone bosons in the external legs. 

4.2 Tadpole terms 

We start with the terms linear in the Higgs fields which will lead to renormalisation of the tad- 
poles. With the tree-level condition on the tadpoles T^o = Txo = 0, field normalisation if it were 
performed does not contribute, we therefore have 

VuneaAo = (57^0 0? + JT^o </)") , (4.7) 



with 






vi Z 5^ + g 

2 

Z" „ I 71 r2 2 \"^'l I / ^2 I ^"Z' 



? + ^C2, + il/|„4) ^ + (-™?, + ^^..)^.^ , (4.8) 

(5T0O ^^00 

= ^(wi ^^ i'2,mi <-^ 7/12) . (4.9) 

The minimum condition requires the one-loop tadpole contribution generated by one-loop dia- 
grams, T°o°^ is cancelled by the tadpole counterterm. This imposes 

6T^o = -r^°P. (4.10) 

T°o°^ is calculated from the one- loop tadpole amplitude for iJ°, T^°^ and /i'^, T^°^ by simply 
moving to the physical basis 



4.3 Mass counterterms in the Higgs sector 

We now move to the mass counterterms induced by shifts in the Lagrangian parameters. We 
need to consider all terms bi-linear in the fields. From the bare matrices M^o, ^1± a-nd M?o 
(Eqs. (|2.6p , (|2.22p - (|2.24p ). we find the corresponding counterterms in matrix form in the basis 
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2 ''2/3 t(i 
^^ia _2 _2 <5tfl r___2 1 _ r7i,r2 -'^Al-i _ Stri 



5m: 



5ML 



6ml + i(44 - l)5Ml, - Ml^slf/if 8m\^ - \s2fiml, - ^s^p'-^ 



42 - \s2p6Mi. '-¥'^p'4f ^^l + 5(44 - i)'5^^lo + ^4o^^,^ 



It is then straightforward to move to the physical fields through the rotation matrices Via) and 
C/(/3), to find the mass counterterms SM^o,SAlfj^,SM^o,SMfjo for, respectively, the pseudoscalar 
Higgs, A^, the charged Higgs H^, and the two CP-even Higgses h^, iJ". A mass mixing between 
these two Higgses, SAlfjof^o is also induced 



6Mlo = slSml + cfiSml - S2p5ml^ - -cjpSMJo + — ^s^^C2/3y^ , 

SMfj± = 5M\o + 5M^± , 

SMJja = c^Sml + s\5m\ + S2a5m\2 

I f , 2 2 2 2 \ \ c 9 ^^7" f \ St0 

+ ^ I 4(C„C^ + SaSp ~ CaSaCpSp) - 1 I M/^o ^520 I 2c2qS20 + S2aC2/3 I — , 

1/^22 22 \ \ 2 ^^7" f \ 5^0 

+ ^ I 4(C„S^ + S„C^ + CaSaCpSp) - 1 1 (51/^0 H ^^2^ I 2c2aS2/3 + S2qC2/3 I — , 

^^'^H°h° = C2Q(5m^2 + -S2q((^TO2 ^ <5m^) 

1 / \ 2 ^^0 / \Stp ^ ^ 

- -^l2s2aC2p + S2pC2ajoMzo-\ ^S2/3 I 2s2qS2/3 - QaQ/J I — ■ (4-12) 

A mass term seems to be induced for the Goldstone bosons as well as a mixing between the 
Goldstones and the corresponding CP-odd Higgs 

SM^o = clSml + slSml + S20Sml2 + -clpSMzo - -MzoslpC2p-— , (4.13) 

Z Z tp 

5Ml± == 5Mlo, (4.14) 

SMqoao = C2pSml2+cpSp{Sml~ Smj)- -C2pS2p6Mzo + Mzoslpcpsp-— , (4.15) 

Z tp 

SMl^^i = SM^oAo ~ M^^cpsp^ . (4.16) 

tp 

It is much more transparent to re-express these mass counterterms by trading-off 6mi^2 and 6mi2 
with our input parameters STa,o , SAl'^o , Stp through 

Sml = c%{s% + 1)^ - 44^ + s%m% - \c2pml. + \slf,{M\. + A'/|o)^ , 

^m?2 = ^-2,(4^ +c|^ - ml. - C2,Mlo^) , (4.17) 
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In terms of ST^a , SM^q , Stp , the mass counterterms of Eq. ()4.12|) and Eq. ()4.16|) write 

SM^o = 6Mq± = -{Ca-pSTHO - Sa-pST^o) , 

SM^OjiO = -{Sa-pSTHO +Ca-pSTho) - 32/3—^-— , 
V Zip 

^^^G±H± = -{Sa-pSTija + Ca-p5Ti^o) ~ S2P f~ ^ ' 



(5M|± = 5M%+5Ml,± 



SMlo = 



--(ca-psl_pSTH» +Sa-pil + cl^p)STho\ +c\_p5M\o +sl+pSMlo 



71 r2 ^^0 

+ S2pS2(a+P)J^J-Z«-r- 1 
t/3 



5Ml 



% = - ( Ca-p{l + sl_p)dTHO + Sa~pcl^_p5Tha j + 4-/3'5M|o + 4+/3<5M|o 



71 r2 '^^/3 
-S2pS2{a+p)Mzo—— , 
t/3 



= -sI-P^Tho + -cl_p5ThO + -S2{a-p)SM'^o - -zS2(a+p)SM'^o ^ f M'^oC2(a-P) 

+ A/|oC2(„+/3))^. (4.18) 



SMfjoho 



It is very satisfying to see that 6Mqo = SM^± is accounted for totally by the tadpole counterterms. 

4.4 Field renormalisation 

We can now introduce field renormalisation at the level of the physical fields without the need to 
first go through field renormalisation in the basis 4>i2i'Pi2^ '/'f 2- ^"^ most generality we can write, 
as in Eq. (|46l) 

ro \ / r'O \ / 7I/2 7I/2 \ / r-o 

^Vo " ^" Uv - Iv ^i^co z)i:^ ) \ ^" 

Lf 1 _'7('-^ 1 — 1 ^|^± ^r±Ff± 1 I "-^ 



TjO \ / TjO \ f 7I/2 7I/2 \ / ttQ 



h' ), '\h' )"\ Zliil, Zli^ \ h- 



hORO ^/i" 



(4.19) 



It is always possible to move basis through Eq. (|4.6p . Field renormalisation will help get rid of 
mixing between physical fields when these are on-shcU and set the residue to 1. 
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4.5 Self-energies in the Higgs sector 

Collecting the contribution of all the counterterms, including wave function renormalisation, the 
renormaliscd self-energies write as 

^G»Aoiq^) = SGO^o(q2) + SM^o^o - ^^SZcoAo - ^{q^ ~ M\,)5Za«G0 
^A"A»iq^) = S^o^o(q2) + SMlo - {q^ - MI,)5Za» 

^G±G±iq^) = SG±G±(g') + SM^^ - qHZG± 

^G±H±[q^) = ^G^H^iq^) + 5Ml^„^ - W5Zg±h± -\{q^~ Mfj^)6ZH±Gi 

^H±H±{q^) = ^H±H±{q^) + SMfji - {q^ - Mfj±)SZH± 

±HOHo{q'-) = SffOHo(g2) + SMfjo - {q^ - MJjo)5Zho 

^HOho{q^) = Si/o^o(g2) + SMfjo^o -^{q^- Mfj„)SZHah« - 5(9' - Ml„)6Zh.H0 

S;,o^o(g2) = E^o^o(g2) + 5Mlo - [q^ - Ml„)SZ^,o 

Note that as we stressed all along, since we are only interested in having finite S'-matrix transitions 
and not finite Green's functions there is no need trying to make all two-point functions finite. For 
instance the diagonal Goldstone self-energies T,QOQo{q'^) and Sg±g±(g2) do not need any field 
renormalisation. Therefore we can set for example 6Zqo = 5Zq± = for simplicity. SZaog° is 
also not needed as it is only involved in the transition of Golsdtone boson to the pseudo-scalar 
Higgs. 

4.6 A°Z° and H^W^ transitions 

The (massive) gauge bosons and the pseudo-scalar mix. This originates from the same part of the 
gauge Lagrangian where the gauge bosons, at tree-level, mix with the Goldstone bosons as in the 
Standard Model, see for example [9]. The latter is eliminated through the usual 't Hooft gauge 
fixing. To wit, from 



I 
9 



J^o = ^i{vid^(b^+v2d'^(b^)W+ + h.c. 



2cw 



(t;i9V; + «2a'V^)^°|o, (4.20) 



we end up with 



l/„ . <5M, 



2 



CGV^^GV ^ _i^sZa±+SZ^i+^^]{iM^^d'^G-W+ + h.c.) 
- ^ (sZgo + SZzozo + ^^ ) Mzo d^G^Zl 



1 



^^z^o^M^oa^G^^ 






\ UzG±ti± + S2/3^) {iMy,±d^H-W^ -1- /i.e.) 



^(.ry St 

- dZoOAO + S2/3 — 
2 V tf3 



- ( SZgoao + S2p^ )Mzod^A''zl . (4.21) 



For the sake of completeness, we have also kept in Eq. (j4.2ip the wave-function renormalisation 
constants of the gauge bosons, namely SZ\y±, SZzozo and SZzo^ (for the Z° — > 7 transition), see 
[S] . The conditions on the latter are the same as in the Standard Model, details are found in [5] . 
The novelty however is that now we have A" — Z^ and H^ — W^ transitions whose self-energies 
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write: 

E^o^o(g2) = j:^,^,^q^) + Mll(sZGOAO+S20^-^Y (4.22) 

SH±v^±(g') = ^H±w±iq^) + ^^(sZG±H±+S2f^Y^y (4.23) 

Note that apart from Stfj the same counterterm SZqojio appears in the G'^A'^ transition. In fact 
there is a Ward identity relating these two transitions. Contrary to what one might see in some 
papers [HI [ITl [18] , the relation is much more complicated for q^ ^ -^lo ^tnd gets more subtle in the 
case of the non-linear gauge. This identity is very important especially that in many approaches 
the transition has been used as a definition for Stfs. The identity can be most easily derived by 
considering the BRST transformation on the ("ghost") operator {0\c^ {x)A^ {y)\0) = 0. Full details 
are given in Appendix [A] We have the constraint 

q^±AOZoiq^)+Mzo±AOGoiq^) = (g^~M|o) ././ f "^'^ a^-^Mlg') (4-24) 






(<?' - M^o) ( 7Jz^z^^!'cHq^) + S2p^ ' 6Zao&> 



J'Q\{<f') and J-liy{q^) arc functions defined in Appendix [Al They vanish in the linear gauge with 
e = 7 = 0. The constraint shows that even in the linear gauge (j^E^02o(g^) + Mzai'A'''G'^>{'f') is 
zero only for q^ = M\o and not for any q^ . We will get back to the exploitation of this constraint 
later. A similar constraint relates also Yifj±-^±{q^) and Y,Q±fj±(q^) 

q'tH^W^{q')+Mw^±H^G^{q') = (g^ - M^±) \/'f '^^ g§^^(g^) 

I \ (47rj^ s^^ tp 

g^HVvil'^) and gPf'^^iq^) are defined in Eq. ^EM), see Appendix |H 

4.7 Renormalisation conditions 

4.7.1 Pole masses, residues and mixing 

Masses are defined as pole masses from the propagator. Moreover this propagator must have 
residue 1 at the pole mass. In the case of particle mixing, the mixing must vanish at the pole 
mass of any physical particle, i.e. at the pole mass. In general in the case of mixing this requires 
solving a system of an inverse propagator matrix with solutions given by the pole masses. For a 
2-particle mixing one has to deal with the determinant of a 2 x 2 matrix which is a quadratic form 
in the self-energies whose solutions are the corrected masses. The equation reads 

= 0. (4.25) 



q - ^4«.trcc - ^h»h"iq )][q - A^HO,troc - ^HOH"iq ) - ^h"H»iq ) 



M/jO tree refers to the tree- level mass. This equations simplifies considerably at one-loop since one 
only has to keep the linear term, or first order in the loop expansion, in the equation. In principle 
the argument that appears in the self-energy two-point functions is the pole mass which might get 
a correction from its value at tree-level. To get the corrections one can proceed through iteration, 
starting from the tree-level masses as argument of the two-point function. Higher order terms in 
the expansion will appear as higher orders in the loop expansion and we do not count them as being 
part of the one-loop correction. A genuine one-loop results for the pole mass, Mj^noop, starting 
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from a tree-level mass -Mi^trcc with Tm^q^) the diagonal rcnormaliscd self-energy is therefore the 
solution of 

q^-Ml^^^-Retu{q^)^0 at g' - Af ^^^^p , (4.26) 

which in the one-loop approximation means 

A^'uoop = A^'trcc + Re%^ (Mltroo) = Mltreo + SMf^ + R^^^^ (A^'trec) ■ (4-27) 

The latter condition will constrain the Lagrangian parameters with 5Mf^ a gauge invariant quan- 
tity. Likewise, at one-loop, the requirement of a residue equal to one, for the diagonal propagator 
and vanishing mixing when the physical particle is on-shell leads to 

i?eS^(M5,,J = with ^^^ = S^(g2)^ 

Ret[^ {Ml,,,,) = Ret[^ {Ml,,,) = t ^ j. (4.28) 

In our renormalisation programme, Eqs. (|4.28p set the field renormalisation constants and avoid 
having to include corrections on the external legs. The field renormalisation constants are therefore 
not necessarily gauge invariant nor gauge parameter independent. 

4.7.2 Renormalisation conditions and corrections on the mass parameters 

As we have explained earlier one needs to fix the counterterms for SM'^a and Stp once tadpole 
renormalisation has been carried through to arrive at finite and gauge invariant S'-matrix elements. 
Taking M^o as an input parameter means that its mass is fixed the same at all orders, we therefore 
set 

SMlo = -ReY^AOAo {M^o) . (4.29) 

Finding a condition to fix Stp is an arduous task that has been debated for sometime. We will 
study many schemes for Stp in Sectional 

The charged Higgs mass is independent of t^, it gets a finite correction at one- loop once M^o is 
used as an input parameter 

Mfi±sioop = ^4±,trcc + Re^„±H±{Mfji^,,„) - Re^AOAo{M%) - i?en^±(M^±) , (4.30) 

we have used 6M^^ = Reli^±{M^±) where n^±(q^) is the transverse 2-point function of the 
W^ following the same implementation as performed in [5] . The finiteness of the corrected charged 
Higgs mass is the first non trivial check on the code as concerns the Higgs sector. 
The sum rule involving the CP-even Higgs masses Eq. (|2.14p is also independent of tp. This sum 
rule gets corrected at one-loop 

A^'oaioop + A4o,iioop = Mlo + M|o + i?eS^o„o(A/2„) + ReE ho h" {M],„) 

+ 7777 ( Ca-pSTno - Sa-pSTffi I - i?eE^o^o (Af^o) - iJeH^o^o (A/|o) . 

(4.31) 

Here also we have used (5Af|o = i^eH^o^o (A^2° ) "^h^re H^o^o {o^) is the transverse 2-point function 
of the Z^ boson, see [3] • Otherwise to predict Af^o ^^ or Mfjo ^^ one needs a prescription on 
dtff, see Eq. (|4.18p . Obviously fixing one of these masses, for instance M^o in particular in analogy 
with M^o, is a scheme for tjs. In this scheme therefore ReT, ^o h° {Mfja) = which sets a gauge 
invariant counterterm for t^, see Eq. (|5.13p . 

14 



4.8 Constraining the field renormafisation constants 

We have introduced through the field renormafisation matrices Zp.Zc.Zs a total of 12 such 
constants, see Eq. (|4.19p . However as argued repeatedly, some of these constants are only involved 
in the transition involving an external Goldstone bosons, i.e. in situations that do not correspond 
to a physical process. Therefore we can give the constants 5Zqo^5Zq±,5Zaoqo,5Z[j±q± any 
value, S'-matrix elements will not depend on these constants. It is therefore easiest to set these 
4 constants to in actual calculations and give them arbitrary values in preliminary tests of a 
calculation of a physical process. 

For the transitions involving physical Higgs particles we just go along the general lines described 
in Section [4. 7. 11 in order to avoid loop corrections on the external legs. In the following, in order 
to avoid too much clutter the masses that will appear as argument are the tree-level masses (or 
the input mass for il/^o). The conditions read 

(4.32) 
(4.33) 
(4.34) 
(4.35) 
i?eEffOfto(M,^o) = 0. (4.36) 

From these we immediately derive 6 out of the 8 field renormafisation constants in the Higgs sector 

5Zao = ReY!j^oAo{M%), (4.37) 

5Zh± = i?eS];,±^±(il4±), (4.38) 

5Zho = ReY!jjoHo{MJjo), (4.39) 

5Zho = ReT:^oAMl-), (4-40) 



i?eS^o^o(Af2o) = 


= 0, 


teS;,±^±(i\4±) = 


= 0, 


i?eS;,„^o(M^o) = 


= 0, 


Rety^oAMlo) = 


= 0, 


RetH-ho{Mlo) = 


= 0, 



SZnoH^ = 2 --^-;;^-^"; „--Ho.o ^ ^^^^^^ 



J.,j2 _ J.J2 



^Zho.o =. 2 - - - '-"^ -"'^" . (4.42) 

When considering a process with A^ as an external legj, in principle it involves the A^ -^ A^ 
transition but also the ^° — > Z'^ and the A" — *■ G^ transitions. The field renormafisation constant 
dZj^o, see Eq. (|4.37[) allows to set the A'^ -^ A^ transition to and moves its effect to a vertex 
counterterm correction. One therefore would be tempted by setting S^ozo(M^o) = together 
with S^O(3o(M^o) = as is done almost everywhere in the literature. In our case this would mean 
that the remaining constant SZqoao could be derived equivalently from one of these conditions. 
However the Ward identity we derived in Eq. (j4.24p imposes a very important constraint. It shows 
that in a general non-linear gauge we can not impose both S^o^o (iW^o) = and S^oqo (M^o) = 0. 
It looks at first sight that this requires that one introduces loop corrections on the external legs 
when considering for example processes with the pseudoscalar Higgs as an external leg. In the 
linear gauge on the other hand this is possible since ^Q^iq^) = 0, we could then adjust SZqoao 
and SZ^OQa to have E^0 2o(M^o) — and S^O(3o(M^o) = 0. Note however that contrary to what 
we encounter in some publications, see for example JTBt \Tf\ . q'^T,AOz"{q'^) + Mz«'S,a°g"{q'^) does 
not vanish for any value of q^ but only for q^ = M^oLJ. 

Let us show how despite the constraint in Eq. (|4.24p we can still avoid one-loop corrections and 
counterterms in the external legs associated with an external pseudoscalar A'^. Of concern to 



*The argument with the charged Higgs is exactly the same, therefore we will not make explicit the detailed 
derivation of the field renormalisation constant &Zq±^± but only quote the result. 

tThe charged counterpart of this identity is also not valid for any q^ as is assumed sometimes, see |18| . 
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us are the transition A'^ — Z^ and A^ — G^ . The idea is that although we can not make both 
YiAozo{M\o) = and EyiO(3o(Af^o) = 0, we wiU try to make the combined contribution to the 
external leg vanish. This combined contribution is pictured in Fig. [TJ 




G" 




+ 




= M 



^xt. leg 



Figure 1: The combined contribution of the A" — Z'^ and A^ — G° transitions 



To the tree-level coupling of the A" to some vertex V, at one-loop the transition A" — G° involves 
the coupling of the tree-level neutral Goldstone to this vertex, Vg while the Z'^ transition involves 
the corresponding vertex V^. The total contribution of Fig. [T] for A^ with momentum q on-shell 
with q^ = M^o writes 



■'^'cxt. leg 



^AOGoiMlo)VG + q.Vz^AOZo{Mlo) 



M% 



Af2 



Vg 



M% 






I]^OGo(Af^o) 4- Mz"^A"Z»iMj^o) 



(4.43) 



In the second step of Eq. (I4.43|l we used another identity that can be readily derived at tree-level 
from the invariance of the Lagrangian under gauge transformationqj. Therefore in order not to 
deal with any correction on the external pseudo-scalar leg we require 



Saogo(MIo) + Mzoi:AozoiMlo) = 0. 



(4.44) 



For this requirement Eq. (j4.44p . which is a renormalisation condition, to be consistent with the 
Ward identity in Eq. (|4.24|) leads to 



S^0 2o(M^o) 



il/. 



z" 



-S^OGo(M^o) — 



1 e^Mzo 



(47r)2 



s^i^J^'SAiMio) 



^2W 



(4.45) 



In particular with J-q'^{M'^o) = in the linear gauge, we can make TiAoz"{M'^a) — I]^0(3o(Af 
0. This condition readily gives 



A<>) 



5Z, 



GO A" — -S2I3- 



,^A«Z<>\^-^A<>) 
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A'L 



zo 



(47r) 



2 „2 
^2W 



S2;3^Ga(^-^Ao)' 



(4.46) 



Since SZaog" only enters in off-shell processes, A'^ off-shell or an external Goldstone boson, there 
is no need to constrain it through some other renormalisation condition. Our aim, as stressed 
repeatedly, is not to renormalise all Green's functions, but only S-matrix elements without the need 
for external leg corrections. The Ward identities that we derived in this section were, numerically, 
checked extensively in our code for various values of q^ including q"^ = M\o and q^ = M|o and 
for different values of the non-linear gauge parameters. Moreover it is thanks to the !Fq]^{M^o) 
contribution in SZqoao that we are able to obtain finite and gauge invariant results for processes 
involving A^ as an external particle. For 5Zg±h± ^ similar derivation gives 



SZg±h± = -S2/3- 



6t^ 






(M^±) 



Mw± 



(47r) 



2 ,— «'(A^^±] 

^2W 



(4.47) 



-f Consider the part of the Lagrangian with the Z^ and the neutral Goldstone G". Before gauge- fixing this 
Lagrangian is invariant under the transformation Zj] —> Z^ -\- 'idi_i,ui, G'' — > G" -I- M^qlo. If the Z" (vector) current 
is V^ and the Goldstone current Vq, that is we have the interaction Z^ .Vz -\- G'^Vq, invariance of the Lagrangian 
implies —idaVz "I" ^zo^G = 0. In Eq. H4.43| l. this implies q.Vz = MzoVq where q is the Z^ momentum. 



16 



With SZqoj^o (and 5ZQ±fj±) all our field rcnormalisation constants are set and defined. 

5 Definitions of tj3 and the tp schemes 

5.1 Dabelstein-Chankowski-Pokorski-Rosiek Scheme (DCPR) 

This scheme, which we will refer to as the DCPR scheme, has been quite popular and is based 
on an OS rcnormalisation scheme in the Higgs sector |14[ I15| working in the usual linear gauge. 
The definition of tp however is difficult to reconcile with an On-ShcU quantity that represents 
a direct interpretation in terms of a physical observable. One first introduces a wave function 
rcnormalisation constant, SZfji, for each Higgs doublet Hi, i.e. before rotation 

H,^{l + ^dZH,)H, 1 = 1,2. (5.1) 

To make contact with our approach and parameters, as concerns wave function rcnormalisation, 
we refer to Appcndix[Bl The vacuum expectation values arc also shifted such that the countcrtcrm 
for each Vi writes 

v^^vAl-^ + -SZh, , (5.2) 

giving 

Stp 5vi Sv2 1, , , , 

-— = -{^ZH^-6ZH^). (5.3) 

tp Vi V2 2 

The DCPR scheme takes ^ = ^ such that in effect 

^-^^\{5Zh._-5Zh,). (5.4) 

tjj is defined by requiring that the (renormalised) A^Z'^ transition vanish at q^ = -Ai^o, therefore 
from 

RetAozo{M\o)=G, (5.5) 

with 

S^02o(g2) = S]^0 2o(g2) + ^.,^p{5Zh, - 5Zh, + 2^) , (5.6) 

4 i/3 

one obtains that 

r, DCPR ^ 

-^ =-- i?eS]^02o(Aflo). (5.7) 

tfi MzoS2l3 

This definition is clearly not directly related to an observable. Moreover 5tp is expressed in terms 
of wave function rcnormalisation constants, sec Eq. (|5.4p . 



5.2 DR Scheme (DR) 

In this scheme the countcrtcrm for tp is taken to be a pure divergence proportional to the ultraviolet 
(UV) factor in dimensional reduction. Guv 

Cuv = 2/(4 -n)--iE + ln(4^) , (5.8) 
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where n is the dinicnsionahty of space-time. In this scheme the finite part of the counterterm is 
therefore set zero: 

-^ =0. (5.9) 

The divergent part can be related to a few quantities not necessarily directly related to an observ- 
able. In the vein of the DCPR approach within the linear gauge where St^ is defined in Eq. (|5.4p . 
solving for SZ^^ — SZ^i leads to the HHW prescription of HoUik, Heinemeyer and Weiglein [TH] , 
see also Eq. (|B.15p . 

r. DR-HHW -, 

-^ =^—{ReY:,^oAMl.)~ReT:HOHo{Ml,)r. (5.10) 

The superscript °° means that only the infinite Cjjv part in dimensional reduction is taken into 
account. A more satisfactory DR scheme can be based on a physical observable. Pierce and 
Papadopoulos [20] have defined Stp by relating it to the divergent part of M'jja — M'^a- Note 
that the sum Mfjo + M^o does not depend on t^ as can be seen from the tree-level sum rule in 
Eq. (PH)) . Hence, see also Eq. (|OT|) . 



5t^^«-^^ 1 /I 



Q(Ca-;3(l + 2sI_p)STho + So_/3(l + 2cl_f,)Sno) 



ip 2s2/3S2(a+/3)A/^o 

+ Re^HOHo{Mfjo) - ReJ^hOhoiMf^o) + C2{a+[})Re^AOAo{MJ^o) - C2{a+l3)ReUzozoiMzo) 

(5.11) 

5.3 An On-Shell Scheme (OSa/jj) with Mho as an input 

In this scheme one takes Mj^o , the largest of the two scalar Higgs masses, as an input parameter. 
This trade-off is operative in the Higgs sector independently of any process. Therefore M^o is no 
longer a prediction but is extracted from a measurement together with Mj^o . As such it does not 
receive a correction at any loop order, Stfj is defined from the constraint 

Rei:H«Ho{M]io)=0, (5.12) 

which leads to 

5i,OS-« _ 1 (.^2 ^2^2 /^*? ,.,2 ^2^2 ST,o 

tp S2/3S2{a-p)M^o \ '^ '^ Vi ' ' V2 

+ ReJ:HOHo{Mfjo)-sl_pRe^AOA"{Mlo)-cl^f,ReUlozo{Mlo)Y (5.13) 

This scheme has been advocated in [T71 [12] and is one of the scheme implemented in SloopS. At 
tree-level, i^ is extracted from the relation defined in Eq. (|2.16p 



MloAffo 



^h = ^ " J J""' "" - (5.14) 



In our numerical examples the input parameters are such that the requirement c^a < 1 is always 
met. In fact given a set M^o, Mzo we generate M^a through a given value of tp. The value M^jq 
is taken as the physical mass at all loop orders, in particular at one-loop it does not receive a 
correction. As pointed out in Section 2, in general with a set A///o , Af^o , Af^o Cja < 1 is not 
guaranteed. With this important proviso, we extract tan/3 (with tan/? > 1) as 



ip 



Mao Mzo + Mho J Ml, + M|o - Mjj, 

V _ (5.15) 

^ Mao Mzo - Mho^M%+MI,-MJj, 
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That this choice might lead to large corrections and large uncertainty can already be guessed by 
considering the uncertainty on tan (3 given an uncertainty on Mfja , M^o , Mza with respectively 
5Mfja , 8Mjyd , 8Mza ■ For clarity let us take SMzo = as would be fit from an experimental point 
of view since M^o is known with an excellent precision from the LEP measurements. We find 

tp ^ Ml„ - M%, Mjj, - MlX M% M% M%, Mfjo M^^ T ' ' 

With typical input parameters in the decoupling limit A/^o ^ Af^o with M^o/Af^o ~ 1 a large 
uncertainty ensues, to wit 



6t0 1 / 5Ml, 5Ml„ 



tp M^o/A/lo-lV A/2„ A/^o 



(5.17) 



Therefore although 5tfj is manifestly gauge invariant one should expect large uncertainty from 
loop corrections. This scheme is similar to the one considered in [6] based on Eq. (|5.14p . 

5.4 At-t- as an input parameter (OS^^^) 

(3 which appears in the Higgs sector relies on the assumption of a basis, only quantities which are 
basis independent are physical quantities [51 [^. The Higgs potential of the MSSM appears as a 
general two- Higgs doublet model if one restricts oneself solely to the Higgs sector. The degeneracy 
is lifted when defining the Yukawa Higgs coupling to fermions. This picks up a specific direction. 
One should therefore define tan/3 from the Higgs couplings to fermions. Since Af^o is used as an 
input parameter assuming one has had access to the pseudoscalar Higgs, it looks natural to take 
a coupling A^ f f . Since couplings to quarks are subject to large QCD radiative corrections the 
best choice is to consider the Arr coupling which is the largest coupling to leptons. 

Ca^ =i spT^zTA ^i^rr-. tp tj^t A with vi ^ vcp ■ (5.18) 

Vi ZM^f± 

This coupling can be extracted from the measurement of the width T a^^ with TOt- the mass of the 
T. Note also that 5Ta^^ = ^dtp/tp so that contrary to the On-Shell scheme based on Mho, OSmhj 
this scheme should therefore not introduce additional large uncertainties assuming of course that 
this decay can be large and be measured precisely. This scheme appears therefore very natural, 
however it has not been used in practice because one has considered it as being a process dependent 
definition set outside the purely Higgs sector which moreover implies that fixing the counterterm 
involves a three-point function. This last argument is unjustified, take for example the G^ scheme 
in the SM where muon decay is used as a trade-off for Af^F* taking advantage of the fact that G^ 
has been for a long time so much better measured than M^r± . The G^ scheme involves four-point 
functions. We find that technically this scheme is not more difficult to implement than a scheme 
based on two-point functions. The full counterterm to Arr involves the G" -^ A'^ shift, the A" 
and r wave function renormalisation constants among other things, we get 

.A.. fSmr 6e cl. SM^^ 1 5M%. 1 1 ~ 

CT ynir e^ 2s\, Af^± 2s\ Ml, ^ 2 ^^ 2tp "" ^ 

-^5Za..^ - f^^^^#^-T:^^M,o^§l(A^|o). (5.19) 

Itfj tp Mza l + t^j Ztt 

Srur, the r mass counterterm, 6e the electromagnetic coupling counterterm, SM\^±zo the gauge 
bosons mass counterterms and the r wave function renormalisation constant SZJ^ ^ counterterms 
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5Ca^^ = CX^\^5^Y + -^j with 



are defined on-shell exactly as in the SM [S] . The fuU one- loop virtual corrections consist of the 
vertex corrections, 5^^^ which contributes a one-loop vertex correction to the decay rate as: 

^pVcrtcx ^ 2ro (5^" . (5.20) 

The latter are made UV- finite by the addition of the counterterm in Eq. (|5.19[) . These virtual QED 
corrections, both vertex and counterterm (from 5mT and (5Z£ ^) include genuine QED corrections 
through photon exchange which are infrared divergent. In our case the infrared divergence can 
be trivially regularised through the introduction of a small fictitious mass. A, for the photon. 
As known, the fictitious mass dependence is cancelled when photon bremmstrahhmg is added. 
Taking into account the latter may depend on the experimental set-up that often requires cuts 
on the additional photon kinematical variables. Therefore it is much more appropriate to take 
as an observable a quantity devoid of such cuts, knowing that hard/soft radiation can be easily 
added. Fortunately for a neutral decay such as this one which is of an Abelian nature, the virtual 
QED correction constitutes a gauge-invariant subset that can be trivially calculated separately. 
The virtual QED corrections to the decay width A^ -^ t^t^ are known [22], they contribute a 
one-loop correction 

^^QED ^ 2ro<5«^^ with 



Li.(i^Ulni±^lni±4-iln^^ + ^''^' 



I + /32 

''"^'l + [3) '"' 2/3 "'!-/? 4'" 1-/3 ' 3 

U2{x) = - / -ln(l-i). (5.23) 

This QED correction only depends on M^o,e,77i,- as it should and does not involve any other 
(MSSM) parameter. Subtracting this QED correction from the full one-loop virtual correction 
in Eq. ()5.19p will give the genuine SUSY non QED contribution that does not depend on any 
fictitious photon mass nor any experimental cut. Our scheme is to require that Stjs is such that 
this contribution vanishes and that therefore A^ -^ t^t~ is only subject to QED corrections. This 
gives 



^^^ - h^--+5^Y -5^'''']. (5.24) 



This definition is independent of the fictitious mass of the photon A used as a regulator. We have 
checked this explicitly within SloopS. 



6 Set-up of the automatic calculation of the cross sections 

All the steps necessary for the renormalisation of the Higgs sector as presented here together with 
a complete definition of the MSSM have been implemented in SloopS. As we will discuss in a 
forthcoming publication [23| the other sectors have also been implemented and results relying on 
the complete renormalisation of the MSSM have been given in [T^ . Since even the calculation of 
a single two-point function in the MSSM requires the calculation of a hundred of diagrams, some 
automatisation is unavoidable. Even in the SM, one-loop calculations of 2 ^ 2 processes involve 
hundreds of diagrams and a hand calculation is practically impracticable. Efficient automatic 
codes for any generic 2^2 process, that have now been exploited for many 2^3 O [25] and 
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even some 2^4 [521 HZ] processes, are almost unavoidable for such calculations. For the clec- 
troweak theory these are the GRACE-loop [5] code and the bundle of packages based on FeynArts 
P5] . FormCalc [55] and LoopTools ^\, that we will refer to as FFL for short. 
With its much larger particle content, far greater number of parameters and more complex struc- 
ture, the need for an automatic code at one-loop for the Minimal Supcrsymmctric Standard Model 
is even more of a must. A few parts that arc needed for such a code have been developed based 
on an extension of |31j but, as far as we know, no complete code exists or is, at least publicly, 
available. Grace-susy [321 is now also being developed at one-loop and many results exist |13j . 
One of the main difficulties that has to be tackled is the implementation of the model file, since 
this requires that one enters the thousands of vertices that define the Feynman rules. On the 
theory side a proper rcnormalisation scheme needs to be set up, which then means extending 
many of these rules to include countcrtcrms. When this is done one can just use, or hope to 
use, the machinery developed for the SM, in particular the symbolic manipulation part and most 
importantly the loop integral routines including tensor reduction algorithms or any other efficient 
set of basis integrals. 

SloopS combines LANHEP [33j (originally part of the package COMPHEP [31]) with the FFL bundle but 
with an extended and adapted LoopTools [llj . LANHEP is a very powerful routine that automati- 
cally generates all the sets of Feynman rules of a given model, the latter being defined in a simple 
and compact format very similar to the canonical coordinate representation. Use of multiplets and 
the supcrpotcntial is built-in to minimize human error. The ghost Lagrangian is derived directly 
from the BRST transformations. The LANHEP module also allows to shift fields and parameters 
and thus generates counterterms most efficiently. Understandably the LANHEP output file must be 
in the format of the model file of the code it is interfaced with. In the case of FeynArts both the 
generic (Lorentz structure) and classes (particle content) files had to be given. Moreover, because 
we use a non-linear gauge fixing condition [9], see below, the FeynArts default generic file had to 
be extended. 



7 tp scheme dependence of physical observables, gauge in- 
var iance: A numerical investigation 

In this first investigation we will restrict ourselves to Higgs observables. Other observables in- 
volving other supersymmetric particles require that we first expose and detail our rcnormalisation 
procedure of the chargino/neutralino and the sfermion sector. This will be presented in [25 . We 
have however presented some results on the tan (3 scheme dependence of a few cross sections that 
are needed for the calculation of the relic density in the MSSM [T5] . 

7.1 Parameters 

To make contact with the analysis of [6] and also allow comparisons we will consider the 3 sets of 
benchmarks points for the Higgs based on [35j . The 3 sets of parameters called mhmax, large fi and 
nomix are as in |35| except that we set a common tri-linear Aj to all sfermions for convenience. 
For each set there are two values of tp, tp = 3, 50. 



7.2 Gauge independence and the finite part of tjs 

If tfj is defined as a physical parameter then St/s must be gauge invariant and gauge parameter 
independent. Our non- linear gauge fixing allows us to check the gauge parameter independence 
of 6tfj and hence tp. Even when two schemes are gauge parameter independent the values of St/^ 
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Parameter 


Value 


Parameter 


Value 


Constant 


Value 


Sw 


0.48076 


■nif. 


0.1057 


nis 


0.2 


e 


0.31345 


rrir 


1.777 


rrit 


174.3 


9s 


1.238 


rriu 


0.046 


rrib 


3 


Mzo 


91.1884 


rrid 


0.046 


Mao 


500 


me 


0.000511 


rric 


1.42 


tp 


3;50 



mhmax 


Value 


nomix 


Value 


large /.i 


Value 


A* 


-200 


M 


-200 


A* 


1000 


A/2 


200 


M2 


200 


A/2 


400 


M3 


800 


M3 


800 


M3 


200 


Mp^ 


1000 


Mp^ 


1000 


Mf, 


400 




1000 
2000+^/i/3 


Af 


1000 


Af 


400 



Table 1: The set of SM and MS SM parameters for the benchmark points. All mass parameters are 
in GeV. We take Mi according to the so-called gaugino mass unification with Mi ™ ^ 



Sci 



are not expected to be the same. It is therefore also interesting to inquire how much two schemes 
differ from each other. Naturally since 6tf} is not ultraviolet finite we split this contribution into 
a finite part and infinite part, the latter being regularised in dimensional reduction, such that 



5tR = 5tr 



Sf^Cuv 



(7.1) 



0. When calculating observables in this scheme we 



The DR schemes have by definition St^''^ 

will also need to specify a scale JI which we associate with the scale introduced by dimensional 

reduction. For the latter our default value isjl = M^o. Our set of non-linear gauge parameters is 

defined as nigs = (a, /?, 5, (D, p, k, e, 7). 

The usual linear gauge, nigs = 0, corresponds to all these parameters set to 0. For the gauge 

parameter independence we will compare the results of the linear gauge to a non-linear gauge 

where all the non-linear gauge parameters have been set to 10, referring to this as nigs = 10. 

To make the point about the gauge parameter dependence it is enough to consider only one of the 
benchmarks points. 



St^ 


nigs = 


nigs = 10 


stT 


nigs = 


nigs = 10 


DCPR 


-3.19x10-^ 


-1.04x10-1 


DCPR 


-0.10 


-0.27 


OSa/„ 


-3.19x10^2 


-3.19x10-2 


OSa/„ 


+0.92 


+0.92 


OSa^^ 


-3.19x10-2 


-3.19x10-2 


OSa.. 


-0.10 


-0.10 


DR-HHW 


-3.19x10-2 


+5.32x10-2 


DR-HHW 








DR-PP 


-3.19x10-2 


-3.19x10-2 


DR-PP 









Table 2: Gauge dependence of Stp at the scale 11 = Ma" for the set mhmax at tp ~ 3. 

As expected we see from Table [2] that only the schemes based on a physical definition of tp are 
gauge parameter independent. Therefore neither DCPR nor a DR manifestation of it based on 
[l9] are gauge independent. Within the physical definitions note that although the divergent part 
is, as expected, the same for all the schemes in all gauges, the finite parts are quite different from 
each other, in particular the OSmh scheme introduces a "correction" of about 30% to t/3. This is 
just an indication that this scheme might induce large corrections on observables. However one 
needs to be cautious, in the same way that the Cuv part cancels in observables, a large finite 
correction could, in principle, also be absorbed when we consider a physical process. Our rather 
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extensive analysis will show that this is. after all, not the case. Schemes where the finite part 
of Stf} is large do, generally, induce large corrections. It is important to note that for the linear 
gauge all schemes give the same Cuv part. Having made the point about the gauge parameter 
dependence, we will now work purely in the linear gauge since some of the schemes introduced in 
the literature are acceptable only within the linear gauge. Therefore in this case the results for 
DR-HHW and DR-PP are the same and will be denoted as DR in what follows. 



7.3 



5tJ'" 



tf3 = 3 


mhmax 


large fi 


nomix 


i/3 = 50 


mhmax 


large fi 


nomix 


DCPR 


-0.10 


-0.06 


-0.08 


DCPR 


+3.42 


+14.57 


+0.48 


OS^ffj 


+0.92 


-1.31 


+0.64 


OSmh 


-385.53 


-2010.84 


-290.18 


OSa^^ 


-0.10 


-0.06 


-0.08 


OSa^^ 


+0.12 


-4.72 


+0.16 


DR 











DR 












Table 3: StL"^ for the Higgs benchmark points. 



First of all let us mention that our numerical results concerning the DCPR and DR schemes agree 
quite well with those of [6] concerning the shifts in tp and the lightest CP-even Higgs mass. Our 
results for OSa/„ follow sensibly the same trend as the scheme defined as the Higgs mass scheme in 
[6]. We see that for small tp DCPR and OS^i^^ give sensibly the same result with a finite relative 
shift of a few percent. For larger i^ the difference is much larger, we notice that OSyi^^ gives much 
smaller shifts. On the other hand the OSj\/„ gives huge corrections for tp = 50 well above 100%. 
As we will see this will have an impact on the radiative corrections on some observables based on 
this scheme. 



7.4 Higgs masses and their scheme dependence 



t/3 = 3 


mhmax 


large /i 


nomix 


t;3 = 50 


mhmax 


large /i 


nomix 


Aq^ = 72.51 








M^o^ = 91.11 








DCPR 


134.28 


97.57 


112.26 


DCPR 


144.50 


35.88 


124.80 


OSj\/h 


140.25 


86.68 


117.37 


OSj\/h 


143.76 


13.21 


124.16 


OSa.. 


134.25 


97.59 


112.27 


OSa.. 


144.50 


35.73 


124.80 


DR /I = Mao 


134.87 


98.10 


112.86 


DR 71 = Mao 


144.50 


35.77 


124.80 


DR 71 = Mt 


134.47 


97.55 


112.38 


DR 71 = Mt 


144.50 


35.77 


124.80 



Table 4: Mass of the lightest CP-even Higgs at one loop in different schemes All masses are in 
GeV. 



We start with the one-loop correction to the lightest CP-even Higgs. Of course, this has now been 
calculated beyond one-loop as the one-loop correction is large, however a study of the scheme 
dependence is important. Moreover this study represents a direct application of the code that can 
be compared to results in the literature. We note that all schemes apart from OSmh s-'^s ™- very 
good agreement with each other for both values of tp. Leaving aside the case of tp = 50 in the 
large n scenario, despite the very large shifts we observed in St J" for the OSmh scheme, the tp 
dependence is much suppressed such that the OSa/^ scheme compares favourably with the other 
schemes. In the case of the correction to the heaviest CP-even Higgs at one loop, by definition 
there is no correction in the OSmh scheme, the other schemes agree with each other at a very 
high level of precision. Moreover especially at high tp the correction is very small. 
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i/3 = 3 


mhmax 


large jU 


nomix 


t0 = 50 


mhmax 


large fi 


nomix 


M]jt = 503.05 








7\/j^^ = 500.01 








DCPR 

OSmh 

OSa.. 

DR 71 = Mao 

DR 71 = Mt 


504.68 
503.05 
504.68 
504.52 
504.63 


501.05 
503.05 
501.05 
500.95 
501.05 


504.21 
503.05 
504.21 
504.08 
504.19 


DCPR 

OSmh 

DR 71 = Mao 
BRJI= Mt 


499.80 
500.01 
499.80 
499.80 
499.80 


498.90 
500.01 
498.91 
498.91 
498.91 


499.85 
500.01 
499.85 
500.01 
499.85 



Table 5: Mass of the heaviest CP-even Higgs at one loop in different schemes. All masses are in 



GeV. The one-loop result is based on the relation Mf^„ = M'^o ^^^^ + ReTj(M^o tree)- 



The mass of the charged Higgs docs not depend on Stp 
pendent, with the counterterm (5A/^± = 6M^± 



therefore the correction is scheme inde- 



SMlo- 



7.5 Higgs decays to SM particles and their scheme dependence 

7.5.1 yl° -^ T^T^, the non QED one-loop corrections 



tp^3 


mhmax 


large /i 


nomix 


r^-^ = 9.40 X 10-3 








DCPR 

OS Mh 

OSa^^ 

DR 71 = Mao 

DR 77 = Mt 


+3.56x10-5 
+6.41x10-3 


+6.51x10--* 
+2.30x10--* 


-8.71x10-'^ 
-7.82x10-3 


+3.94x10--* 
-2.66x10-5 


-7.37x10"" 
+4.56x10-3 


+5.18x10--* 
+9.67x10-5 


tp = 50 


mhmax 


large /i 


nomix 


r'^'^ = 2.61 X 10" 








DCPR 

OSmh 

DR 71 = Mao 
DR 71 = Mt 


+3.45x10-* 
-4.03x10* 


-1.21x10-2 
-3.00x10-2 


+2.01x10" 

-2.09x102 



+4.92x10-* 

+4.75x10-* 


+3.35x10-2 
-3.03x10* 


-1.66x10-2 
-3.44x10-2 



Table 6: Corrections to the decay A" 
All widths in GeV. 



at one loop without the universal QED correction. 



We now study the non QED corrections to the decay width A" -^ t^t~, see Section \U^ for our 
benchmark points. By definition there is no correction in the OS^^^ scheme. Many interesting 
and important conclusions can be drawn from Table [Sj First of all we note that the scheme 
dependence is quite large here. After all this is an observable which is directly proportional to St p. 
In fact the difference between schemes can be accounted for by 2(5t^ read off from Table [3l For 
this decay, the OSi\ifj scheme is totally unsuitable, for i^ = 3 the correction are of order 100%, 
whereas for i^ = 50 the one-loop correction is an order of magnitude, at least, larger than the 
tree-level. Especially for i^ = 3 in DR the scale dependence is not negligible. For example with 
Jl = mt in DR the correction is of order ^ 1% and 5% for Ji = Mao. The corrections are much 
smaller in DCPR being at the per-mil level. The scale dependence is much smaller for tp = 30 
and the corrections in DR are now smaller than in DCPR. Note also that in the large ^ scenario 
the corrections arc large. 
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7.5.2 H'^ -^ T^T , the non QED one-loop corrections 



</3 = 3 


mhmax 


large fi 


nomix 


r^ -^ = 9.35 X 10-3 








DCPR 


-1.09x10-4 


-7.96x10-5 


-1.09x10-4 


OSmh 


+6.28x10-3 


-7.91x10-3 


+4.47x10-3 


OSa.. 


-1.45x10-'' 


-7.09x10-5 


-1.01x10-4 


DR 71 = Mao 


+5.08x10-4 


+3.24x10-4 


+4.17x10-4 


DR 77 = Mt 


+8.57x10-5 


-9.75x10-5 


-4.52x10-6 


ti3 = 50 


mhmax 


large ^ 


nomix 


V^ = 2.61 X 10" 








DCPR 


+3.54x10-^ 


+2.02x10" 


+4.31x10-2 


OSa/„ 


-4.03x10^ 


-2.09x102 


-3.03x101 


os^.. 


+9.52x10-3 


+ 1.94x10-3 


+9.55x10-3 


DR 71 = Mao 


-2.59x10-3 


+4.94x10-1 


-7.00x10-3 


DR 77 = Mt 


-2.04x10-2 


+4.76x10-1 


-2.49x10-2 



Table 7: Corrections to the decay _ff" 
All widths in GeV. 



at one loop without the universal QED correction. 



Similar conclusions can be drawn from the study of the non-QED corrections to i/" -^ t^t~ , see 
Table [T] The QED corrections for this decay can be implemented as in [22] . The only difference 
is that now there is a correction also in the case of the OSa^^ scheme. But as expected this 
correction is very small for both values of tp. Note that for t^ = 50 the DCPR scheme gives 
very large corrections in the large n scenario. For this process we have not taken into account 
the one-loop correction to Mfjo , since as we have seen this correction is very small for all schemes 
and also because one is much too far from the tt threshold, Mfjo ^ 500 GeV ^ 2tot-, where this 
effect can play a role. 



7.5.3 i/o ^ Z°Z° and A" ^ Z°h° 



i/3 = 3 


mhmax 


large fi 


nomix 


F'^ ^ = 8.97 X 10-3 








DCPR 

OS Mh 

OSa^^ 

DR 71 = Mao 

DR 71 = Mt 


+ 1.59x10-2 
+ 1.40x10-2 
+ 1.59x10-2 
+ 1.57x10-2 
+ 1.58x10-2 


-6.32x10-3 
-4.00x10-3 
-6.32x10-3 
-6.44x10-3 
-6.32x10-3 


+8.47x10-3 
+7.12x10-3 
+8.47x10-3 
+8.32x10-3 
+8.44x10-3 


ti3 = 50 


mhmax 


large fi 


nomix 


T^^ = 6.40 X 10-5 








DCPR 

OSmh 

os^„ 

DR 7t = Mao 
DR p = Mt 


+2.18x10-5 
+1.01x10-2 
+3.02x10-5 
+3.05x10-5 
+3.09x10-5 


-5.14x10-4 
+4.66x10-3 
-4.65x10-4 
-4.77x10-4 
-4.76x10-4 


+3.89x10-5 
+7.81x10-4 
+3.97x10-5 
+4.01x10-5 
+4.05x10-5 



Table 8: Corrections to the decay H'^ -^ t^t at one loop. All widths in CeV. 

i/" -^ Z^Z^ was studied by [5^ where a large correction was found. Wc confirm here, see Table|51 
that a large correction is indeed induced with the one-loop result of the same order if not exceeding 
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both at tfj = 3 and tp ~ 50 the tree- level result. This larger correction is not due to the scheme 
dependence since in this process the latter is very small whereas one sees a large correction with 
all the schemes. The correction is large because the benchmark points with M^o = 500 GeV are 



O7O7O 



in the decoupling regime where if -^ Z 'Z" practically vanishes at tree-level. The H"Z^Z 
proportional to c^_q ~ M^o/M^o, the coupling is therefore almost induced at one loop without 
the 1/il/^o suppression. Here again because Alffo ^ 2AIzo the one-loop correction on M^o is 
negligible. Very similar results and conclusions can be drawn for the process A'^ -^ Z'^h^, see 
Table El 



tfi = i 


mhmax 


large fi 


nomix 


V'^'^ = 9.03 X 10-3 








DCPR 

OSmh 

DR 71 = Mao 
DR 71 = Mt 


+2.42x10-2 
+2.23x10-2 
+2.50x10-2 
+2.48x10-2 
+2.41x10-2 


+3.86x10-3 
+6.20x10-3 
+3.86x10-3 
+3.74x10-3 
3.87x10-3 


+ 1.68x10-2 
+ 1.55x10-2 
+ 1.64x10-2 
+ 1.67x10-2 
+ 1.68x10-2 


ti3 = 50 


mhmax 


large ^ 


nomix 


r' ^ = 6.30 X 10-5 








DCPR 

OSa/„ 

DR 71 = Mao 
DR 71 = Mt 


+2.39x10-5 
+1.00x10-3 
+3.48x10-5 
+3.51x10-5 
+3.30x10-5 


+8.75x10-'' 
+5.97x10-3 
+9.26x10-4 
+9.12x10-4 
+9.12x10-4 


+4.31x10-5 
+7.74x10-4 
+4.39x10-5 
+4.43x10-5 
+4.47x10-5 



Table 9: Corrections to the decay A'^ -^ Z'^h^ at one loop. All widths in GeV. 



8 Conclusions 

The use of the non-linear gauge has allowed us, for the first time, to quantitatively and quali- 
tatively, study different proposals for the ubiquitous parameter tan/3 and its effect on the Higgs 
observables, both the physical Higgs masses as well as their decays. Our first preliminary con- 
clusion is that the scheme based on the extraction and definition of tan (3 from a decay such as 
A^ -^ T^T^ is by far the most satisfactory. Not only is this definition directly related to a physical 
observable and therefore is gauge independent, the functional dependence of the physical width in 
tan/3 is linear and is the same independently of the value of the pseudoscalar Higgs mass. More- 
over the definition is clean once we subtract the universal gauge invariant QED correction. The 
scheme is also most pleasing and satisfactory since it is the one where the observables we have 
studied show the least corrections, leading therefore to a stable prediction. On this last count the 
DR scheme performs almost just as well. However the widely used DR scheme extracted from 
the A^Z^ transition is not gauge invariant and therefore terribly unsatisfactory from a theoretical 
point of view. In the non-linear gauge with a general gauge-fixing set of parameters, the parameter 
gauge dependence shows up already at one-loop, whereas it has been known that the scheme fails 
even in the linear gauge but at two-loop [7]. A gauge independent DR scheme such as the one 
proposed in [^^ is most satisfactory. A scheme based on the usage of Mjjo as an independent 
parameter from the Higgs sector leads to too large corrections in most of the observables we con- 
sidered so far. We therefore propose that the decay A° -^ t^t~ be used as a definition of tan/3. 
This choice assumes that this decay will one day be measured with high enough precision but this 
depends much on the spectrum of the MSSM. Were it not for the unambiguous extraction of the 
full QED corrections, the decay of the charged Higgs to tv may also qualify as a suitable input 
parameter, see j36j for prospects on the measurement of this decay. Apart from the discussion on 
gauge invariance and the issue of the scheme dependence for tan /3, we have shown how a complete 
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one-loop renormalisation of the MSSM can be automatised and have given resuhs and details as 
concerns the Higgs sector which is the first step in a successful implementation of this programme. 
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Appendices 

A The Ward-Slavnov- Taylor identity for the transitions A^Z^ 

and A^G^ 

There is an identity relating the A'^Z'^ and A°G° transition. This is most useful for q^ = M^o- 
Contrary to what one might see in some papers, the relation is much more complicated for q^ ^ 
Af^o and gets more subtle in the case of the non-linear gauge. 

The identity can be most easily derived by considering the BRST transformation on the ( "ghost" ) 
operator {0\c^{x)A'^{y)\0) = 0. We find 

<5BRs(0|c^(x)A°(y)|0) = (0|(<5brsc^(x))AO(2/)|0) - {0\c^{x)iSBRsA'{y))\0) =0, (A.l) 

with 

SbrsA° = -^{c+H-+c-H+) + ^c''{c^^fih" + s^^pH°), (A.2) 

Z S2W 

and Sbrsc^ = B^ . (A.3) 

Therefore. 



{Q\B^{x)A\ym + |(^(0|c^(x)c+(y)7?-(y)|0) + (0|c^(x)c-(y)if+(y)|0) 

c„_^(0|c^(x-)c^(j/)/i°(2/)|0) + s„_^(0|c^(x)c^(j/)i7"(2/)|0)') - . (A.4) 



521^ 



At tree-level, there is no vertex involving c^c^H^. Using the equation of motion of the B field, 
we obtain a relation for the following Green's functions (external legs are not amputated): 

d^{0\Z^{x)A°{y)\0) + Mzo{0\G%x)A\y)\0) 



S2W 

e 



S2W 

In a diagrammatic form, we have 
1 1 



e{0\h'\x)G"ix)A''{y)\0)+j{0\H''{x)G"ix)A^y)\0) 
c,-^(0|c^(a;)c^(2/)/iO(z/)|0) +.s„_;j(0|c^(a;)c^(y)iJ°(2/)|0) 1=0. (A.5) 



i e 



q2 - Ml, g2 _ ^/2 



iqf, X Z^' --^O--^ A° + Mzo X G° --^O— ^ A^ 



q^ - M\a S2W 



Ux 0^o°--^^° + 7X Ogo-^^°) 



and obtain the relation 



' Ca-fi X C^ --^0^,0 +Sa_/3 X C^ --^0^0 ) , (A.6) 



_. „ ^ -Z .r^C^ I „ _ w -Z 

^2 - Mia S2W 



q^Y.AOz^'{q^) + ilVE^OGo(q2) ^ _(g,2 _ A/|„)^^ ( gx O^o'"^ ^° + 7^ Ogo"^ ^° ) 

S2W V / 

+ {q^ - Mlo)-^(cc,-,3 X c^ -^Ofo +so,^0 X c^ -^O^^'o ) . (A.7) 

S2W V / 

With the following vertices 

£ D _£i^,2//..+/3/.°-C„+^i/°V°G"' (A.8) 

S2W V / 

C^' D ^^f(.„_0-6>"-(c„_,+7)if"')c^c^, (A.9) 

S2W \ J 
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we calculate all the "lollipops" 



-i^M^S2f}So,+pBo{q^M^o,M^zo) , (A.IO) 

S2W 
iSM^S20C^+pBoiq^Mfj„,A4a) , (A.ll) 

eMzo 
I 

-^(3ho = ~i^^^{c^-p+l)Bo{q^Mfjo,Mlo), (A.13) 

5214^ 





-^^0 




-^AO 


c^- 


-0^: 



(s„_/3 ~ ?)So(g', Af^o, M|„) , (A.12) 



with 

Boiq^MlMl) = Cuv- I dx ln{A{q^ , M^ Ml)) , (A.14) 

Jo 
A{q^,M^,Ml) = q^x^ - (g^ + Af| - Af2)a; + A/S . (a.15) 

Wc finally obtain the identity 

q^J^AOZoiq') + MzoJ:AOGo{q') = tAs ^"^^ f ('^^ " M|o)s2;3^^1(g2) + {q^ - M%)Tl^Hq^] 

with T~^\{q^) = jCa,+f}Boiq\ Af^o, M|o) - «„+/jBo(g2, M^o, Af|o) , 
^f'^(g2) = &„_/3Bo('z', Af^„, Af|„) + 7s„_^So(g2, Af^o, Af|o) 

+ ^S2(„_0) (i?o(g', Af^o, A/|o) - Bo{q\Ml..Mlo)) . (A.16) 

To implement this formula into SloopS and check it numerically, we need to introduce the tadpole 
part in FormCalc and we define S'^"^ the self-energy without tadpole: 

g^E^XoVC?') + MzoT.'X^a«{<l^) + Mzo6T = ^^^^((g^ _ m|o)s2;3^M + (g^ - M%,):f!^) , 
where ST = ——{sa-pSTno + Ca-i}5Tu°) ■ (A. 17) 

S2wMzo 

We remark on some simplifications in the functions JF for specific choices of the non linear gauge 
parameters 

^^'l(e = 0,7 = 0) = 0, (A.18) 

J^aAie = c^+p,l = s^^,) ^ -.2(.+,)y^ '^"^"U(g^M|o,A^|o)j' ^^'^^^ 

J^'J{e = s^^fi,^ = -Co,^p) = 0, (A.20) 

In terms of renormalised self energies, 

E^02o(g2) = j:^-i^^q^) + M^(^sZaoA^.+s20^-^), (A.22) 

E^OGo(g2) = E^^o^Go(g2) + 5A/l„c„ - ^q'SZaoAo - ^(g' - Aflo)5Z^OGO , (A.23) 
with (Eq. I4l8)) 



(JA^loGO = 5r--S2/3Mjo-^, (A.24) 
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we obtain the following constraint on the renormalised two-point functions 



q^±AOZoiq^) + Mzo±AOG"{q^) - (g^-A/|o) .. ., ^ ^% 2/3.FSl(g') (A.25) 






(9' - M%')i7J^7^^'c^il^) + ^2p^- 5Zag) 



Note that in this identity 5T and more importantly 5Zqoa° drop out. 

The derivation of the identity for the charged Higgses follows along the same steps. We only quote 
the result 

2 V(47rj-^S2^ tp 

with the functions Q^w il^) ^^^ Scf'^iQ-"^) defined as 

ShW (9^) = S{s2pSa+p - C^CQ_/3)Bo(g^, M^±, A^ho) " '^(c2/3Sa+/3 - S^Sq_^)Bo (?^ , M^± , M|fo ) 

+pc^Bo(g',M^±,M2o), 

-pclrBo{q^,M^±,Mlo). (A.26) 

B Wave function renormalisation constants before rotation 

In our approach field renormalisation was performed on the physical fields, or better said, after 
rotation to the /i", H^, A°, G", H^, G^ basis. We could have applied field renormalisation on the 
components of the doublets Hi,H2, Eq. (j2.3p . To make contact with some of the early papers 
[in [T51 [TO] on the renormalisation of the Higgs sector we therefore introduce the most general 
field renormalisation on the components of iJi , i/2 ■ We define 

^'i\ - ( ^ ^Zi \(^'^\ mi) 



*" ^0 " ZS.. Z^i \ <^t 



(B.2) 



:ix^(IS)(i)- 

As explained in the text, these constants are immediately transformed into the set of matrices 
Zp.ZctZs- Or we can go from the set Zp,Zc,Zs to the set defined by Eqs. (jB.lHB.Bp . For 
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example, 



SZgo = CfjSZ^o + Sp5Z^o + CpSp{5Z^o^o^ + 5Z^a^o) 
SZgoj^o = C[3S0{5Z^o - SZ^o) + cfjSZ^o^a - s'p6Z^a^o 
SZjiOQO = CpSi3{5Z^o - SZ^o) + c|(5Z^o^o - s'pSZ^o^o 
dZjiO = s'^pSZ^o + C^f^SZ^o - Cf3S[3{SZ^o^o + SZ^o^o) 

5Zg± = cl5Z^± + sl6Z^± + C[3Sp(6Z^±^± + 5Z^±^±) 
5Zg±h± = cpsi:s{5Z^± - SZ^±) + cl5Z^±^± - slSZ^±^± 
6Zh±g± = cpSfi{SZ^± - SZ^±) + c}SZ^±^± - s}SZ^±^± 
SZh± = slSZ^± + clSZ^± - cpsp{SZ^±^± + SZ^±^±) 

SZffO = C^SZ^o + sldZ^o + CaSa{SZ^o^a + SZ^o^o) 
^ZnOhO = CaSa{SZ^o - 5Z^o) + C^6Z^o^o ~ slSZ^o^o 
SZh«H« = CaSa{5Z^o - 5Z^a) + C^^SZ ^0^0 - s\5Z^a^Q 
5Zh« = s\5Z^o + c\5Z^o - CaSaiSZ^o^o + SZ^o^o) 

SZnOhO + SZjiOHO = (Cq, - S^){SZ^o^o + SZ^o^o) + 2CaSa{SZ^o - SZ^o) 

5ZH0h° — SZfiOHO = SZj,Oj,0 — iZ^OAfi 



(B.4) 

(B.5) 

(B.6) 

(B.7) 



Our rcnormalisation conditions in Eq. (j4.36p on S^,; will turn into 

i?eI]^o^o(A/|o) = s^(5Z^o +c^(5Z^^o - c/3S^((5Z^o^o +^^^0^0), (B.8) 

i?eE^±^± (M|r± ) == s%8Z^± + c^(5Z^± - c^s^(5Z^±^± + Z^±^±) , (B.9) 

-ReE^o^o(A/^n) = c^fJZ^o + s^^Z^o + CqSq(5Z^o^o + SZ^o^o) , (B.IO) 

■ReSlo^o(M^o) = slSZ^o + clSZ^o - c^SaiSZ^o^o + SZ^o^o) . (B.ll) 



In fact in [HI [HI [19] only two renormalisation constants are introduced, one for each doublet 
through 

H,^{l + hzH,)H, z = l,2. (B.12) 
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This means that 



6Z^o = 6Z^o = SZ^± = 6Zh, , 
SZ^o^o = SZ^o^o ^SZ^±^± =0 i^j. (B.13) 

Since wave function rcnormalisation is applied on the doublets it also contributes a shift to Vi. 
Another shift on this parameter is also applied, Vi — > Vi~Svi as to all other Lagrangian parameters. 
Compared to our shift 5vi, we have 

Sv,=Sv,--SZHiVi. (B.14) 

Note that with only SZh, and SZh.,, m view of Eqs. (|RT0)) - ((RTT|) and Eq. (|RT3)) we have 

SZh, - 5Zh, = - ^ (ReKoho {Ml, ) - i?eS;,„^o (M^o )") . (B.15) 



31 



References 



[1] Y. Okada, M. Yamaguchi and T. Yanagida, Prog. Theor. Phys. 85 (1991) 1; 
ibid. Phys. Lett. B262 (1991) 54; 

J.R. Ellis, G. Ridolfi and F. Zwirner, Phys. Lett. B25 (1991) 83; 
ibid. Phys. Lett. B262 (1991) 477; 

H.E. Habcr and R. Hcmpfling, Phys. Rev. Lett. 66 (1991) 1815; 
A. Yamada, Phys. Lett. B263 (1991) 233. 

There is a long list of papers on the one-loop corrections to the lightest Higgs mass. It is fair 
to add to the list the seminal paper, written with the assumption that the top is not heavier 
than 60 GeV: 
S.P. Li and M. Sher, Phys. Lett. B140 (1984) 339. 



[2] 
[3] 
[4] 

[5] 

[6] 
[7] 
[8] 

[9] 

[lo; 
[11 

[12 

[13 

[14 
[15 
[16 
[17 
[18 
[19 

[2o; 

[21 
[22 
[23 
[24 



A. Djouadi, Phys. Kept. 459 (2008) 1, hep-ph/0503173 



S.P. Martin, hep-ph/9709356 



For a review on the corrections at the two-loop order, see for example 



S.P. Martin, Phys. Rev. D71 (2005) 016012, hep-ph/0405022 



S.P. Martin, Phys. Rev. D75 (2007) 055005, hep-ph/0701051 



R.V. Harlander, P. Kant, L. Mihaila and M. Steinhauser, Phys. Rev. Lett. 100 (2008) 191602, 
arXiv:0803.0672 [hep-ph]. 



A. Freitas and D. Stockinger, Phys. Rev. D66 (2002) 095014, |hep-ph/0205281 

Y. Yamada, Phys. Lett. B530 (2002) 174, ,hep-ph /011225f, 

S. Davidson and H.E. Haber, Phys. Rev. D72 (2005) 035004, [Erratmn-ibid. D72 (2005) 



099902], hep-ph/0504050 



G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato and Y. Shimizu, 
Phys. Rep. 430 (2006) 117, hep-ph/0308080. 



F. Boudjema and E. Chopin, Z. Phys. C73 (1996) 85, hep-ph/9507396 



F. Boudjema, A. Semenov and D. Temes, Phys. Rev. D72 (2005) 055024, hep-ph/0507127 



N. Baro, F. Boudjema and A. Semenov, Phys. Lett. B660 (2008) 550. larXiv:0710.1821l [hep- 
ph]. 

J. Fujimoto, T. Ishikawa, Y. Kurihara, M. Jimbo, T. Kon and M. Kuroda, Phys. Rev. D75 
(2007) 113002. 



A. Dabelstein, Z. Phys. C67 (1995) 495, hep-ph/9409375 



P.H. Chankowski, S. Pokorski and J. Rosiek, Nucl. Phys. B423 (1994) 437,|hep-ph/9303309 
A. Dabelstein, Nucl. Phys. B456 (1995) 25,|hep-ph/9503443 



L.Y. Shan, Eur. Phys. J. C12 (2000) 113, hep-ph/9807456 



H.E. Logan and S.F. Su, Phys. Rev. D66 (2002) 035001, 'hep-ph/0203270 



S. Heinemeyer, W. HoUik and G. Weiglein, Phys. Rept. 425 (2006) 265, hep -ph/041221 4 
D. Pierce and A. Papadopoulos, Phys. Rev. D47 (1993) 222, hep-ph/9206257 



F. Boudjema and A. Semenov, Phys. Rev. D66 (2002) 095007, hep-ph/0201219 



M. Drees and K.I. Hikasa, Phys. Lett. B240 (1990) 455, [Erratum-ibid. B262 (1991) 497]. 

N. Baro and F. Boudjema, paper in draft form, to appear soon. 

G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato and Y. Shimizu, 

Phys. Lett. B559 (2003) 252, |hep-ph/0212261, 

G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato and Y. Shimizu, 



Phys. Lett. B559 (2003) 252, |hep-ph/02 12261; 

G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. Shimizu and 



Y. Yasui, Phys. Lett. B571 (2003) 163, hep-ph/0307029 



32 



G. Bclangcr, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, Y. Kurihara. K. Kato and 

Y. Shimizu, Phys. Lett. B576 (2003) 152,'hcp-ph/0309010' 

F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. Kurihara, Y. Shimizu and 



Y. Yasui, Phys. Lett. B600 (2004) 65, |hep-ph/0407065 
F. Boudjema, J. Fujimoto, T. Ishikawa, T 



Kaneko, K. Kato, Y. Kurihara, Y. Shimizu 
S. Yamashita and Y. Yasui, Nucl. lustrum and Meth. A534 (2004) 334, hcp-ph/0404098^ , 

[25] A. Denner, S. Dittmaier, M 



|hep-ph/0307"l93| 

A. Denner, S. Dittmaier, 

hep-ph/030l"l89| 



Roth and M.M. Weber, Phys. Lett. B575 (2003) 290, 
M. Roth and M.M. Weber, Phys. Lett. B560 (2003) 196, 



You Yu, Ma Wen-Gan, Chen Hui, Zhang Ren- You, Sun Yan-Bin and Hou Hong-Sheng, Phys. 

Lett. B571 (2003) 85, |hcp-ph/03 06036; 

Zhang Ren- You, Ma Wen-Gan, Chen Hui, Sun Yan-Bin, Hou Hong-Sheng, Phys. Lett. B578 



[26] 



(2004) 349, |hep-ph/0308203 

F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. Kurihara and Y. Shimizu, 

Nucl. Phys. Proc. Suppl. 135 (2004) 323, hep-ph/0407079; 

F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato, Y. Kurihara, Y. Shimizu and 



Y. Yasui, |hep-ph/0510184 

[27] A 



Denner, S. Dittmaier, M. Roth and L.H. Wieders, Phys. Lett. B612 (2005) 223, 
[hep-ph/0502063 ; 



A. Denner, S. Dittmaier, M. Roth and L.H. Wieders, hep-ph/0505042 



[28] J. Kiiblbeck, M. Bohm, and A. Denner, Comp. Phys. Commun. 60 (1990) 165; 
H. Eck and J. Kiiblbeck, Guide to FeynArts 1.0, Wiirzburg, 1991; 
H. Eck, Guide to FeynArts 2.0, Wiirzburg, 1995; 



T. Hahn, Comp. Phys. Commun. 140 (2001) 418, hep-ph/0012260 



[29] T. Hahn and M. Perez- Victoria, Comp. Phys. Commun. 118 (1999) 153, |hep-ph/9807565[ 
T. Hahn, hcp- ph/0406288, .hep-ph/0506201, 

[30] T. Hahn, LoopTools, http://www.feynarts.de/looptools/. 



[31] T. Hahn and C. Schappacher, Comp. Phys. Commun. 143 (2002) 54, hep-ph/0105349 



[32] J. Fujimoto et al, Comput. Phys. Commun. 153 (2003) 106, hep-ph/0208036 



[33] A. Semenov, hep-ph/9608488 



A. Semenov, Nucl Inst. Meth. and Inst. A393 (1997) 293; 
A. Semenov, Comp. Phys. Commun. 115 (1998) 124; 
A. Semenov,, hep-ph/0208011; 
A. Semenov. larXiv:0805.0555i [hep-oh]. 

[34] [CompHEP Collaboration], E. Boos et al., Nucl. Instrum 



Meth. 



hep-ph/0403113 



A. Pukhov et al., "CompHEP user's manual, v3.3" 

[hep-ph/9908288, 

http: //theory. sinp.msu.ru/comphep/. 



A534 (2004) 250, 
Preprint INP MSU 98-41/542 (1998) 



[35] M.S. Carcna, S. Heinemeycr, C.E.M. Wagner and G. Weiglcin, piep-ph/9912223 

[36] [ATLAS/CMS Collaboration], Martin Flechl et al, ATL-PHYS-CONF-2007-018, 
[arXiv:0710. 1761^ 4 [hcp-ph]. 



33 



